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ABSTRACT 


Under a NASA MSFC SBIR Phase II effort an analysis has been developed for liquid flows with 
vapor bubbles such as those in liquid rocket engine components. The analysis is based on a 
combined Eulerian-Lagrangian technique, in which Eulerian conservation equations are solved for the 
liquid phase, while Lagrangian equations of motion are integrated in computational coordinates for 
the vapor phase. The novel aspect of the Lagrangian analysis developed under this effort is that it 
combines features of the so-called particle distribution approach with those of the so-called particle 
trajectory approach and can, in fact, be considered as a generalization of both of those traditional 
methods. The result of this generalization is a reduction in CPU time and memory requirements. 
Particle time step (stability) limitations have been eliminated by semi-implicit integration of the particle 
equations of motion (and, for certain applications, the particle temperature equation), although 
practical limitations remain in effect for reasons of accuracy. The analysis has been applied to the 
simulation of cavitating flow through a single-bladed section of a labyrinth seal. Models for the 
simulation of bubble formation and growth have been included, as well as models for bubble drag 
and heat transfer. The results indicate that bubble formation is more or less "explosive": for a given 
flow field, the number density of bubble nucleation sites is very sensitive to the vapor properties and 
the surface tension. The bubble motion, on the other hand, is much less sensitive to these 
properties, but is affected strongly by the local pressure gradients in the flow field. In situations where 
either the material properties or the flow field are not known with sufficient accuracy, parametric 
studies can be carried out rapidly to assess the effect of the important variables. Future work will 
include application of the analysis to cavitation in inducer flow fields. 
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OBJECTIVE: 

DEVELOPMENT OF AN ANALYSIS FOR LIQUID FLOWS WITH VAPOR 
BUBBLES (SUCH AS THOSE IN BEARINGS, SEALS, AND PUMPS) FOR 
USE IN DESIGN OF LIQUID ROCKET ENGINE COMPONENTS 


APPROACH: 

USE A COMBINED EU LERIAN-LAGRANGIAN ANALYSIS 

- CONTINUOUS (LIQUID) PHASE TREATED BY SOLVING EULERIAN 
CONSERVATION EQS. (N-S EQS.) WITH DISCRETE PHASE SOURCE TERMS 

- DISCRETE (VAPOR BUBBLE) PHASE TREATED BY INTEGRATING 

LAG RANG IAN EQUATIONS OF MOTION IN COMPUTATIONAL COORDINATES 


LAGRANGIAN ANALYSIS - METHODOLOGY 


• GENERALIZATION OF PARTICLE TRAJECTORY AND PARTICLE 
DISTRIBUTION MODELS 

=* REDUCED CPU TIME AND STORAGE REQUIREMENTS 

• STABLE INTEGRATION OF LAGRANGIAN EQUATIONS OF 
MOTION IN COMPUTATIONAL COORDINATES 


NO PARTICLE TIME STEP LIMITATION OTHER 
THAN FOR ACCURACY 
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LAGRANGIAN ANALYSIS - FEATURES 


LAGRANGIAN DISCRETE PHASE ANALYSIS INTERFACES 
WITH EULERIAN CONTINUOUS PHASE ANALYSIS 
VIA ONE SINGLE SUBROUTINE 
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GENERALIZED ARRAY ADDRESSING AND INDEX / VARIABLE 
NUMBERING 

=> LAGRANGIAN ANALYSIS CAN BE HOOKED UP TO 
A VARIETY OF EULERIAN CODES ("FLOW SOLVERS") 

CODE IS MODULAR W.R.T. PHYSICAL MODELS 


APPLICATION TO VAPOR BUBBLES 


MODELS NEEDED FOR BUBBLE FORMATION, 
GROWTH / COLLAPSE AND MOTION 


CRITICAL PHYSICAL PARAMETERS: 

- SURFACE TENSION 

- VAPOR PRESSURE - TEMPERATURE RELATION 
(CLAUSIUS-CLAPEYRON EQUATION) 
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BUBBLE FORMATION 

• NUCLEATION DUE TO LOCAL P DROP; 

USE HOMOGENEOUS NUCLEATION MODEL 
(KATZ AND BLANDER) 


/ NUCLEI x _ / 2a x 1/2 r 

" (tst) 


-16jta 3 

3kT(P - P L ) : 


NUMBER OF NUCLEATION SITES PER UNIT VOLUME 
AND TIME 

EQUILIBRIUM VAPOR PRESSURE 

PRESSURE IN LIQUID 

TEMPERATURE OF THE LIQUID 

SURFACE TENSION 

BOLTZMANN CONSTANT 

MOLECULAR WEIGHT 

MOLECULAR NUMBER DENSITY IN LIQUID 


BUBBLE GROWTH 


• AFTER THE INITIAL GROWTH PERIOD (WHICH IS INERTIA 
CONTROLLED) THE PRESSURE INSIDE THE BUBBLE IS 
DETERMINED BY THE YOUNG-LAPLACE EQUATION 


£ - P l =2o/R 


• THE BUBBLE TEMPERATURE IS ASSUMED TO BE 

EQUAL TO THE SATURATION TEMPERATURE AT THE PRESSURE IJ 

• BUBBLE GROWTH IS DETERMINED BY THE ENERGY EQUATION 
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ENERGY EQUATION FOR BUBBLE 


q = m CpTp + m h, g 
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Cp : VAPOR-PHASE SPECIFIC HEAT 

h fg : LATENT HEAT OF EVAPORATION 

m : BUBBLE MASS GROWTH RATE ("EVAPORATION RATE") 

T : RATE OF CHANGE OF BUBBLE TEMPERATURE 
P 

* / ^ \ "* 

t p=UL u - v '’ 


HEAT TRANSFER RATE TO BUBBLE 


q s Nu ji D k (t - T p ) 


: NUSSELT NUMBER 
: BUBBLE DIAMETER 

: THERMAL CONDUCTIVITY OF THE LIQUID 
: LIQUID TEMPERATURE 
: BUBBLE TEMPERATURE 
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BUBBLE EQUATION OF MOTION 




E *a_P d 2 X 

F = m 1+ ^p TT 

l 2 P pJ 



P : 

LIQUID DENSITY 



P p = 

VAPOR DENSITY 



V 

ADDED MASS COEFFICIENT (=1 FOR STOKES FLOW) 



m : 

BUBBLE MASS 



X : 

BUBBLE POSITION VECTOR 
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F : 

FORCE ON BUBBLE 

- DRAG FORCE 

- BUOYANCY 

• PRESSURE GRADIENT 

J 
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DRAG COEFFICIENT 

FOR BUBBLES 


C D ( = Tt) =14.9 /Re 0,78 Re >2 

npa z u z 


= 16 /Re Re <2 

GOOD FIT TO NUMERICAL SOLUTION 

COMPARES WELL WITH EXPERIMENTS, WIDELY USED 


FOR SOLID SPHERES 

C D = — [1 +0.15R( 
Re 

= 0.438, Re > 1000 


[1 +0.15 Re 0687 ] Re <1000 


NOTE: FOR SAME a AND u, DRAG ON A BUBBLE 
IS SMALLER (FACTOR 6 FOR Re = 900) 
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COMPUTATIONAL TEST CASE 


SINGLE BLADE LABYRINTH SEAL 

APPROXIMATES ONE BLADE OF THE LOX PREBURNER 
BOOST PUMP IMPELLER SEAL IN THE SSME 

FLOW CONDITIONS: 


INLET PRESSURE : 3.67 MPa (532 psi) 
INLET TEMPERATURE : 130 K (234 °R) 


DENSITY 

PRESSURE DROP 


1088 kg/m 3 (AG7.9\bli\?) 
1.47 MPa (213 psi) 
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COMPUTATIONAL GRID 




V 
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NEAR-WALL RESOLUTION IN THE GAP = 1.0 x 10' 4 (= 0.015 |im) 



COMPUTATIONAL GRID 
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AXIAL VELOCITY 









BUBBLE TRAJECTORIES 
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